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Abstract 



\ In the testable Flipped SU(5) x U(l)x model with TeV-scale vector-like particles from F-theory 

£NJ ' model building dubbed as the J--SU (5) model, we study the vector-like quark contributions to B 

■ physics processes, including the quark mass spectra, Feynman rules, new operators and Wilson 

I coefficients, etc. We focus on the implications of the vector-like quark mass scale on B physics. We 

find that there exists the sbZ interaction at tree level, and the Yukawa interactions are changed. 
Interestingly, different from many previous models, the effects of vector-like quarks on rare B decays 

X ■ 

$_i , such as B — > X s j and B — > X s £ + £ do not decouple in some viable parameter space, especially 

when the vector-like quark masses are comparable to the charged Higgs boson mass. Under the 
constraints from B — > X s ^f and B — > X s £ + £~, the latest measurement for B s — > can be 

explained naturally, and the branching ratio of B s — > £ + £~j can be up to (4 ~ 5) x 10~ 8 . The 
non-decouling effects are much more predictable and thus the J--SU(5) model may be tested in 
the near future experiments. 

PACS numbers: 12.15.-g, 12.15.Lk, 12,15.Ff, 14.20.Mr, 12.39.-x 
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I. INTRODUCTION 



Supersymmetry provides a natural solution to the gauge hierarchy problem in the 
Standard Model (SM). In the supersymmetric SM (SSM) with R-parity under which 
the SM particles are even while the supersymmetric particles (sparticles) are odd, the 
SU(3)c x SU(2) L x U(l)y gauge couplings can be unified around 2 x 10 16 GeV |l|, the 
lightest supersymmetric particle (LSP) such as the neutralino can be a cold dark matter 
candidate |2|, yj, and the electroweak (EW) precision constraints can be evaded, etc. Es- 
pecially, the gauge coupling unification strongly suggests Grand Unified Theories (GUTs). 
However, in the supersymmetric SU(5) models, there exist the doublet-triplet splitting prob- 
lem and dimension-five proton decay problem. Interestingly, these problems can be solved 



elegantly in the Flipped 577(5) x U(l)x models |4|-|6| via missing partner mechanism jg]. 
On the other hand, string theory is the most promising candidate for quantum gravity, and 
it can unify all the fundamental interactions in the Nature. However, the string scale is at 
least one-order larger than the conventional GUT scale. 

To solve the little hierarchy problem between the traditional GUT scale and string scale, 
two of us (TL and DVN) with Jing Jiang have proposed the testable Flipped SU (5) x U (l)x 
models, where the TeV-scale vector-like particles are introduced [3] ■ Such kind of models can 



be constructed from the free fermionic string constructions at the Kac-Moody level one 
and locally from the F-theory model building 
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111 ], and is dubbed as J r -5 , t/(5) [ll|. In 



particular, these models are very interesting from the phenomenological point of view 11]: 
the vector-like particles can be observed at the Large Hadron Collider (LHC), proton decay 
is within the reach of the future Hyper-Kamiokande "L^j and Deep Underground Science 



and Engineering Laboratory (DUSEL) 13( experiments 14|, |15J , the hybrid inflation can be 



naturally realized, the correct cosmic primodial density fluctuations can be generated 16| . 
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With no-scale boundary 



and the lightest CP-even Higgs boson mass can be lifted 
conditions at SU(5) x U(l)x unification scale 19[, two of us (TL and DVN) with James 
Maxin and Joel Walker have described an extraordinarily constrained "golden point" {20 1 
and "golden strip" 21] that satisfied all the latest experimental constraints and has an 
imminently observable proton decay rate lj]]. For a review of the recent progresses, see 



Ref. 



22]. 



Interestingly, the vector-like quarks in the F-SU^) model predict rich phenomenology 
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on low energy processes. If the model is treated seriously, constraints from electroweak 
parameters such as U, S,T and Rb, R c and B processes should be taken into account. We 
also would like to point out that the J r -5'?7(5) model has no Landau pole problem and then 
is very different from the other simple SM extensions in quark sector (also see the next 



Section) |23|, and the 3x3 SM-like quark mixing matrix is now replaced by a 5 x 5 one and 
then is no longer unitary, and there exists the tree-level sbZ interaction, which will play an 
important, even dominant, role in some parameter space for rare B decays. 

Thanks to the efforts of the B factories and LHC, the exploration of quark-flavor mixing 
is now entering a new interesting era. It is well known that the rare B decays induced by 
the flavor changing neutral current (FCNC) only occur at loop level in the SM and then 
are sensitive to new physics. Thus, the rare radiative, leptonic and semi-leptonic B meson 
decays are valuable in testing the SM at loop level and probe new physics. On the theoretical 
side, the rare B inclusive radiative decays B — > X s ^ and B — > X s £ + £~{£ = e,/i) as well as 
the exclusive decays B s — > u + u~ and B s — > £ + £~^ have been studied extensively at the 



leading 
models 



ogarithm 
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order (LO) 24j and high order in the SM 25] and various new physics 



27( | . On the experimental side, B — > X s "y and B — > X s £ + £~(£ = e,n) have 
been measured and the latest upper bound on B s — > is achieved 28|. By comparing 

the predictions with experimental measurements, we will present some constraints on the 
parameter space in the P-SUfi) model. 

The first task of this work will be deriving the quark mass spectra and Feynman rules. We 
stress that the Feynman rules which not be presented in previous studies are used not only in 
B physics but also in research of all low energy processes. B physics constraints on the model 
is the second task of this work, we will concentrate our attention on the vector-like quark 
contributions to B physics, in particular, the contributions from the new operators induced 
by tree-level FCNC. We will show that the sbZ interaction can be generated at tree level, 
and the Yukawa interactions are changed, new operators 0' g and O' 10 in effective Hamiltonian 
should be introduced. We will demonstrate that different from many previous models, the 
effects of vector-like quarks on rare B decays such as B — > X s 7 and B — > X s £ + i~ do not 
decouple in some allowed parameter space, especially when the vector-like quark masses 
are comparable to the charged Higgs boson mass. Within the constraints from B — > X s ^f 
and B — > X s £ + £~, and the latest measurement for B s — > will be explained naturally, 

and the branching ratio of B s — > £ + £~^ can be up to (4 ~ 5) x 10~ 8 . Because the non- 



3 



decouling effects are very predictable, the F-SU (5) model may be tested in the near future 
experiments. 

This paper is organized as follows. We present a brief description for the TeV-scale F- 
SU(5) model and derive all the Feynman rules for our calculations in Section |TT1 We discuss 
the implications of vector-like quarks on B physics in Section IIHI Our numerical results are 
presented in Section IIVt and Section |V] is the summary. 

II. THE F-SU(5) MODEL AROUND THE TEV SCALE 

To achieve the string-scale gauge coupling unification in the F-SU(5) model, we intro- 
duce the vector-like particles which from complete Flipped SU(5) x U(l)x multiplets. The 
quantum numbers for these additional vector-like particles under the SU(5) x U(l)x gauge 
symmetry are [3] 

XF=(10,1) , YF = (TO, -1) , 
Xf =(5,3) , Ff = (5,-3) , 

XI = (1,-5) , Yl = (1,5) . (1) 

To avoid the confusion in the following discussions, we change the convention in Ref. [t| a 
little bit. It is obvious that XF, YF, Xf, Yf, XI, and Yl are standard vector-like particles 
with contents as follows 

XF = (XQ, XD C , XN C ) , YF = (YQ C , YD, YN) , 
Xf = (XU,XL C ) , Yf = (YU c ,YL) , 

XI = XF , Yl = YE C . (2) 

Under the SU(3)c x SU{2)l x U(1)y gauge symmetry, the quantum numbers for the extra 
vector-like particles are 

XQ = (3,2,1) , YQ c = (Z,2,-\) , 
XU=(3,1,^), YU C = (3,1,-f) , 
XD=(3,1,~) , W = (3,l,i) , 
XL = (l,2,-±) , YL C = (l,2,i) , 
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XE = (1,1,-1) , YE° = (1,1,1) , 

xiv = (1,1,0) , yiv c = (1,1,0). (3) 

At the GUT scale the superpotential is given by 

W GUT = Y^F i F j h + Y^Fjrh + Y^lJ j h + ^hh + Y k ^ ( f> k HF j 
+ Y'fXFFjh + Y'^XFfjh + Y v Y"FpCfh + Y'fXlfjh 
+ Y"fljXfh + Y"f?(f) k HXF + Y 2D XFXFh + Y' 2D Y FY Fh 
+ Y 2Uu XFXfh + F ,2l/iy FFF//i + Y 2E XlXfh + F /2E F/F//i 
+ MjFjYF + MffjYf + Mj%y/ 



+ M XFYF + M b XfYf + M b XlYl , (4) 

where z is the generation indices. The first line is the SSM superpotential, the second line is 
the Yukawa mixing terms between the SM fermions and vector-like particles, the third and 
fourth lines are the SM-like superpotential for vector-like multiplets, and the fifth and sixth 
lines are bilinear mass terms. After the SU(5) x U(l)x gauge symmetry breaking down to 
the SM gauge symmetry, we obtain the superpotential as follows 

Wew = (Xg - yJi)UnQj ■ H d + Y$ v U?Qi ■ H n - Y^N^L, ■ H u 

- Y?E%L ■ H d - Y' 3 D (XD C Q ■ H d + D]XQ ■ H d ) + Y^U^XQ ■ H u 

- Y- Uv XN c L ■ H u + Y- Uu XU c Q ■ H u - Y- Uu N^XL ■ H u - Y- E XE C L ■ H d 

- Y] E E]XL ■ H d - 2Y 2D XD C XQ ■ H d - 2Y >2D YDYQ C ■ H u 

+ Y 2Uu XU c XQ ■ H u - Y 2Uu XN c XL ■ H u - Y' 2Uu YUYQ c ■ H d 
+ Y' 2Uu YNYL c ■ H d - Y 2E XE C XL ■ H d - Y >2E YEYL C ■ H u 

- 2M) \d c -YD + Q ■ YQ C + WYN] + M 2 [U C YU + L ■ YL C ] + M]E C YE 



- 2M 4 [XD C YD + XQ ■ YQ C + XN C YN] + M 5 [XU C YU + XL ■ YL C ] 

+ M 6 XE C YE . (5) 

At low energy, the sparticles decouple rapidly when Ms increases. Note that the LHC 
already put strong constraints on squark masses around 1500 GeV, we will concentrate on the 
contributions from new vector-like quark multiplets XU, YU C , XD, and YD C for simplicity. 
At first glance these multiplets seem to be similar to the fourth and fifth generation quarks, 
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but indeed (XU, YU C ) and (XD, YD C ) are vector-like. This makes them very different 
from the fourth and fifth generation quarks. The down-type quark mass matrix is 



M, 



i) 



' (YR + YR 


)v d 


(YR 


+ YR] 


Vd 


(YR + Yi 


)v d 


YRv d 


-2M{ 


(YR + Yg 


)Vd 


(YR 


+ YR] 


v d 


(YR + Y32 


)v d 


YRv d 


-2M\ 


(YR + Y° 


)Vd 


(YR 


+yr: 


v d 


(^33 + ^33 


)v d 


YRv d 


-2M\ 


YRv d 




YRv d 




YRv d 




2Y 2D v d 


-2M 4 


v 2M\ 






2M\ 




2Ml 




2M 4 


-2Y' 2D v 



(6) 



uj 



and the up-type quark matrix is 



Mr, 



1 11 Uu 


1 21 u u 


yUv 

1 31 u u 


Y ,r Rv u 


Ml 


1 12 u u 


1 22 Uu 


1 32 u u 


Y% v v u 


Mf 


Y Uv v 
13 u u 


I 23 V u 


33 u u 


Y% v v u 


Ml 


ynUv , 
1 1 u u 


yiiUv . 


1 3 u u 


Y 2U »v u 


M 5 


-2M\ 


-2M\ 


-2Ml 


-2M 4 


ynUv 



(7) 



where v u and v d are the vacuum expectation values (VEVs) for H u and H d . These two 
matrixes can be diagonalized by unitary matrices U and V, 



V]M D U d = diag.[m d ,m s ,m b ,m dx ,m dy ], 
V^MuU u = diag.[m u ,m c ,m u m Ux ,m Uy }. 



(8) 



Thus, the quark mixings are described by a matrix V = U\jJ d . From Eqs. (jSJ) and ([7j), 
we can see that the mass matrices of the down-type quarks and up-type quarks are related 
to each other, implying that the Yukawa couplings are different from those in the SM. In 
the Feynman gauge the Feynman rules for charged W boson, Goldstone boson, and charged 
Higgs boson with quarks uidjX + (x — W, G, h) and for Z boson djdiZ needed in our 
calculations are given as follows 
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i^=r[gl(l,j)P L + g x R (l,j)PR}, (X = W, Z) 



V2 

■ 9 

V2 



[g x L (l,j)P L + 9 x R (l,j)P R ], (x = G, h) 



(9) 
(10) 



where 



4 



/ 4 U 

m=l 



u: mi u^\ 9n(i,j) = v: 5i v¥, 



(11) 



9 Z L {hh 



9L(h3) 



9%{hh 



9 h L (h3) 



9 h R ^h 



1 



COS 9yy 
1 

y2 cos 9w 



1 - ^sin 2 e w ) s ij -ufW 



d u d 



2 

3 



~ sin 2 6 w 5 ij + v; 5i vf i 



(12) 



E Y k u X ki uT j + 2r^Kr 5i c/^ 



i fc,m=l 
4 



E (^mfc 
, fc,m=l 



m w 

1 iSJ^tc* - ^t^W 



Vd 

mw ' 



(13) 



E y^v: h uT + 2Y' 2 ^; 5 ^ 



/2D-,r*5i TT 5j I 



^fc,m=l 
4 



E (Y° k + Y k l)V; kj U™ -2Y' u »V^Uf 



(14) 



v k,m=l 

Because the vector-like particles do not change U(1)em interaction, the interactions of pho- 
ton and quarks are still the same as those in the SM. From the above mass matrices we can 
see that the TeV-scale J r -S'6 r (5) model has two points for rich physics to be explored: 

• Since the quark mass matrices are not the same as two Higgs doublet model 
(2HDM) [271 or the Minimal Supersymmetric Standard Model (MSSM) [26|, the loop- 
level FCNC will be changed by the Yukawa interactions, and then may change the 
prediction of process b — > s^y significantly. 

• The last terms in Eqs.f llip -( !T4|) . which we call the "tail terms", will cause the tree- 
level FCNC processes induced by b — > s£ + £~ and then the stringent constraints on the 
model parameter space will be expected. 



III. IMPLICATIONS ON B PHYSICS 

Apart from the directly search for the light vector-like quarks at the LHC, another way 
to test the F-SU^) model is to measure their effects on low energy processes such as rare 
B decays. 



A. Effective Hamiltonian 



The starting point for rare B decays B -» X s 7, B ->■ X s £ + £~, B s -> £ + £~ and B s -> £ + £~^ 
is the determination of a low-energy effective Hamiltonian obtained by integrating out the 
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heavy degrees of freedom in the theory. For b — >■ s transition, this can be written as 

/"f 10 

n cS = ~Kv»El c iM°ito + C 'M°M , (is) 

V 1 i=i 



where the effective operators Oj are same as those in the SM defined in Ref. 24j. The 
chirality-flipped operators O- are obtained from Oi by the replacement 75 — > —75 in quark 
current. It is obvious that 0' g 10 can be got directly from the tail terms in the Feynman rules 
of the ,F-iS77(5) model. A few remarks follow on the operators and Wilson coefficients: 

• As mentioned in introduction, the three generation quark mixing matrix is replaced 
by a 5 x 5 matrix U^Ud and then is non-unitary. In our analyses we take a reasonable 
assumption that the deviation from unitary is not large. Otherwise, the tree-level 
FCNC will modify significantly the low energy processes such as Z — > bb and B s — > 

n + n~. 

• Since the Wilson coefficient C^ipiw^) = — y^y* — 1 is always a good approximation 
in J 7 -SU(5) model, and the coefficients of four quark operators Cj(/Xf,) (i — 1,3 — 6) 
depend actually on the value C 2 (mw), the contributions from the four-quark operator 
matrix elements to effective coefficient Cg (//&) can not be ignored and have the same 
expressions as the SM. 

• The coefficient of operator 0' 2 = (sc)v+A{cb)v-A, for example, is proportional to the 
elements of quark mixing matrix or £/J l denoted the mixings between the ordinary 
quarks and vector-like quarks. Thus, it can be reasonably set to be much smaller than 
O(l), and the contributions from the four-quark primed operators to Cg (/^) and 
C g ' e ^(fib) can be neglected safely. This means 

d$ f M = C' 9>10 (m w ) , (16) 

which receive contributions mainly from the tree-level diagrams, loop diagrams for 
b —> 57, and box diagrams. We also neglect the operator 0' 7 contribution. 

• For b — > 57, the new contributions mainly come from the new type Yukawa interac- 
tions, and for b — > s£ + £~, the new contributions mainly arise from the new operators 
O' 

^9,10 ' 



s 



B. Analyses in B Physics Calculations 



In the J r -5't/(5) model the contributions to operators Oi {% = 1 — 10) and 0' g 10 can 
be encoded by the values of the coefficients C, and C- at the matching scale my/. In this 
Section, we will present the Wilson coefficients at the matching scale and decay widths for 
some rare B decays. We keep both new physics contributions and the SM results at the LO 
for consistency. 

• The Wilson coefficient Cj at the matching scale is 

CV = ^E{^0^*(i.%r(i.3)-B(*0^r(i,%2(i,3) 

VtbVts i= i m b 

TT) 

+ g?(i,2)[C(xi)g2(i,3) - -^I>(x,)ffg(<,3)] 

+ ^g h L *(i,2)[C( yi )g h L (i,3) - ^D( yi )g h R (i,3)]}, (17) 

where Xi = Tn^./rn^ and yi = mf l ./mj i+ . For cross check, using the loop functions 
given in the appendix and the CKM matrix unitarity condition, one can easily obtain 
the predication Cj M (mw) = A(xt) + B(x t ) +x t [C(x t ) + D(x t )} which is consistent with 



that in Ref. [24|. Furthermore, Cj receives a large non-decoupling contribution not 
only from top quark as in the SM but also from the up-type vector-like quark loops at 
the electroweak scale. The non-decoupling effects are unique and will be demonstrated 
in next Section. 

The Wilson coefficient Cg at the matching scale is 

P(x t )-Q(x t ) 

C $ = r^F, K4Q(£i) 

sin t) w 

9-jv /7*52/T53 i 

— — Sill U\y ) 



a em VtbVts 4 

5 



+ 77177 IE [R(xi)g¥*(h2)g¥(hZ) + s^)^*(z\%£(z,3) 

VtbVts U=3 

+ E^T(x,) [^*(z,%«(z,3)+^(i,2)^,3) 

i=l m Ui 

+ ^( yi )rfr(i,2)^(t,3)]J + ^. (18) 

Note the first part related to P(x t ) and Q(x t ) from the box diagrams and the effective 
vertex b — > sZ* at loop level have the same expression as those in the SM, while the 
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second part denotes the interaction at tree level enhanced by a large factor The 
last part comes from the effective vertex b — > 57* at loop level for consistency. The 
contribution from one-loop matrix element of the operator O2 is also included as in 



the SM [24J. Moreover, the Wilson coefficients C10, C' Q , and C[ at the matching scale 
are 

P(x t ) - Q(x t ) , 2tt 1 Uf*V? 

sin 2 ^ + a em A V tb V t * s ' UJj 

1 9 t/*52t/53 

C 9 = (^-sin 2 ^ — ^r^, (20) 

, 2tt 1 y; 52 \zj 3 

The contributions from loop diagrams to C 9 10 can be neglected safely 
Branching Ratios 

Considering that the Wilson coefficients do not separate into the SM and new physics 
parts easily and new operators are introduced, we need to list some explicit expressions 
for the branching ratios of B decays as follows 

1. B^X sl 

The inclusive B — > X s ^ rate is the most precise and clean short-distance informa- 
tion that we have, at present, on AB = 1 FCNCs. The new contributions mainly 
come from the new type Yukawa interactions to operator O7. The calculation of 
the branching ratio is usually normalized by the process B — > X c eTT e , so we get 



Br(5 ^ X sl ) = Br ex (B X c eir e ) 1 % g °? J C 7 e// M| 2 . (22) 





2 6a 






2 


t/(z) 



Here z = and /(z) = 1 — 8z 2 + 82 s — z 8 — 24z 4 lnz is the phase-space 
factor in the semi-leptonic B-decay. From the formula of C7 in Eq. (fTT|) and 
the corresponding coefficients in Eqs. (llip - (114jl . we can see that if we sum the 
flavor indices from 1 to 5 in Eqs. (fTTT)-(fl4l). Cj will be exactly the same as the 
five generation 2HDM. In our numerical calculation we will compare both results 
in these two models, since it will show clearly the implications of the new type 
Yukawa interactions in the F-SU (5) model. 
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2. B->x s e+e~ 

Since the new operators 0' g and 0' 1Q contribute to B — > X s £ + £~ and the exclusive 
decays, the analytical expression of invariant dileptonic mass distribution is found 
to be similar to the SM as follows 
dT(B -»■ X s £ + £~) 



ds 



G 2 F m ! n 2 

768tt 5 eml 
x (4|C 7 e// | 2 



v»v u 

1 + 2 -) 



* |2, 



1-s 



(|C 9 e// | 2 +|^| 2 )(l + 2 S ) 



(23) 



+ (ICiol 2 + |C( | 2 )(1 + 2a) + 12Re(C e / f C c 9 ff *)} 

where s = (pe+ + pi-) 2 /m 2 .. Also, we use the normalization process B — > X c eV e 
to get rid of large uncertainties due to m\ and CKM elements as in Eq. ( |22l . 
3. B s — > ji + n~ 

The purely leptonic decays constitute a special case among exclusive transitions. 
It is strongly helicity suppressed and only receives contributions from two axial- 
current operators Oio and O' 10 in the models we studied. The decay width is 
given by 



n, 2 r 2 



I |2 
10 1 ) 



(24) 



r(B s /iV") = IW fk m B s m^\c l0 - c 

where /b s is the decay constant for B s determined by (01^7^7561-85) = —ifB q Pfi- 
The factor k denotes the non-zero width difference of the S s -meson system effect 
on the branching ratio of the B s — > fi + fi~ decay and it reads 31] 

1 + |r Bs ^ Ar Ar s 



K 



(25) 



where AT S is the difference between the decay widths of the light and heavy 
B s mass eigenstates and tb s is the B s mean lifetime. The parameters »4.Ar is 
related to the effective B s — > lifetime t^+^- and depends sensitively on new 

physics. 

4. B.-ti+e-'y 

The exclusive decay can be obtained from the inclusive decay b — > sl + £~^, and 
further, from b — > s£ + £~. To achieve this, for £ = e,/i we just attach photons to 



any external quark lines in the Feynman diagrams of b — > sV 
rate is 



301 ] . The decay 



dT 

ds 



a 



3 l 2 G f 



-v tb v t * s 



m 



B s 



■8(1 



\K\ 



\L\ 



\M\ 2 + \N\' 



(26) 
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where s = p 2 /m 2 Ba is normalized dileptonic mass squared, and 
K 
L 



\- \ [C 9 e/ V) + C' % \G x [j?) - 2C e 7 "(^G 2 (p 2 ) 



m 



P 



M 



N 



ni 



B 3 



m 



B 3 



with Gi and being the form factors 



(27) 



IV. NUMERICAL RESULTS 



Since additional vector like quark introduced in the model, there are many new input 
parameters appear in Wilson coefficients C 7 , C 9 , C w , Cg, C 10 . These parameters are not 
independent and constrained by conditions Eq. ([HD. As the first study on B physics in 
the model, we will not scan the parameter space completely, but focus on the implication 
of mass scale of the vector-like quark on B physics, this will give us the most important 
information of the model. Thus in the numerical study we scan the mass m Ux in the range 
180 GeV ~ 2000 GeV, and m Uy in the range 40 ~ 60 GeV heavier than m Ux . As for other 
parameters, we use the shooting method to randomly generate 5x5 unitary matrix V u and 
U u , then use the CKM matrix to get the V^, Ud to let mass of down-type quark matrix satify 
the Eq. (Ej). Note that to take in account impact of the non-zero width difference of B s 
system [33| on the branching ratio of B s — > we use y s = 0.088 ± 0.014 3l|. We also 

use the following experimental constraints from B physics: 

1. In the model with three generation quarks, the CKM matrix unitarity is already used 
in the calculations of the loop-level FCNC induced rare B decays. Therefore for con- 
sistency, in the model we study the constraints on CKM matrix element measurements 
are not from rare B decays but from tree-level B decays 34J as shown in Table HI 



2. To see the implications o: 



on the rare B decays 



28 



the vector-like quark multiplets, we use the following bounds 



33] 



Br{b ceV e ) = (10.74 ± 0.16) x 10^ 2 
Br(B X s7 ) = (3.06 ± 0.23) x 10~ 4 
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TABLE I: The CKM matrix elements constrained by the tree-level B decays. 





absolute value 


relative error 


direct measurement from 




0.97418 ± 0.00027 


0.028% 


nuclear beta decay 


V us 


0.2255 ±0.0019 


0.84% 


semi-leptonic K-decay 


Vub 


0.00393 ± 0.00036 


9.2% 


semi-leptonic B-decay 


v cd 


0.230 ±0.011 


4.8% 


semi-leptonic D-decay 


Vcb 


0.0412 ±0.0011 


2.7% 


semi-leptonic B-decay 


Vtb 


> 0.74 




(single) top-production 



Br(B -> X s £ + r) = (4.5 ± 1) x 10" 6 , 

Br(B s -> fj+fj,-) < 4.5 x 10" 9 (95%C.L.) . (28) 

3. Other input parameters are the same as those in the SM, except for tan /3 and the 
charged Higgs boson mass m^+. In our numerical calculations we scan the two pa- 
rameters randomly and choose two typical points (tan — 2, rrih+ = 3000 GeV) and 
(tan (3 = 40, m h + = 500 GeV) for the demonstration. 




200 300 400 500 600 700 BOO 9001000 200 300 400 500 600 700 BOO 9001000 

m Ux (GeV) m Ux (GeV) 



FIG. 1: Comparison of B — > X s ^f versus m Ux in the J--SU(5) model (red cross) and 2HDM (green 
triangle) . 
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The numerical results of B — > X s "y as a function of the vector-like quark mass are dis- 
played in Fig. [TJ For the comparison, Fig. [1] also shows the results of the five-generation 
2HDM. From this figure one can see some features clearly: (i) The new physics effects de- 
couple when the charged Higgs boson is very heavy. However, for a much heavier charged 
Higgs, the branching ratio of B — > increases with m Ux in the J--SU(5) model while is 
almost independent on the extra quark mass in 2HDM, indicating the large non-decoupling 
effects; (ii) Unlike the 2HDM where the large tan (3 is preferred if the charged Higgs boson 
mass is at the EW scale, the small tan (3, which is excluded in 2HDM, is still survived in the 
J-'-SU (5) model; (iii) It is clear from the left plot of this figure that the branching ratio can 
be much bigger than the detection result when m Ux getting close to the charged Higgs boson 
mass. So the detection results of B — > X s ^ can give stringent constraints on the F-SU^) 
model. The tendency of the figure can be understood as following: 



• Cj determined by Eq. ( ITT]) in both J r -5'L r (5) model and 2HDM [27] will approach to the 
SM value when the charged Higgs boson is much heavier than EW scale. Nevertheless, 
the contributions from the fourth and fifth generation up-type vector-like quarks in 
2HDM can be suppressed by small V 5t and V Al due to the unitarity condition of 5 x 5 
matrix; 

• Because the summed indices are only from 1 to 4 in the J r -S'L r (5) model, the unitary 
condition of the CKM matrix can not be maintained. When the vector-like particle 
mass approaches to the charged Higgs boson mass, the suppression from 5x5 CKM 
mixing matrix will be released and then the non-decoupling effects will be sizable. In 
fact, the non-decoupling effects are a very special part of the F-SU^) model at EW 
scale and can be tested at the LHC and other B physics detectors. 

Fig. [2] shows the branching ratio of B — > X s £ + £~ versus B — > X s ^ in the J r -5't/(5) model. 
Clearly, both processes will give stringent constraints on our model. Especially, most part 
of the points are excluded when the charged Higgs boson is several hundred GeV, leaving 
a narrow part in the parameter space. Similar phenomenology can be seen in Fig. [3] which 
shows branching ratios of B s — > fi + n~ versus B — > X s £ + £~ . The non-decoupling effects can 
be stringently constrained by the experiments as expected. Here we should emphasize that 



the upper bounds from the Tevatron and the first LHCb constraints 



331 ] . which are about 
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FIG. 2: Branching ratio of B — > X s £ + £ versus B — > X s j in the J--SU(5) model. 



one order of magnitude above the SM expectation, as well as the recent CDF results of 



B s — > detection 35[ can be explained naturally. It is interesting to see that there is an 

approximate linear relation between branching ratios of B s — > and B — > X s £ + £~. In 

fact, we find that in the allowed parameter space with Ud — Vj, the dominant contributions 
to both processes come from C; and C\{i = 9, 10). From Eqs. (1181) to (121j) . we can easily 



draw the conclusion that the branching ratios are nearly proportional to \C 



10 1 



m h . = 3000 GeV I »« 
tanf, = " 





Br(B^X s ri) 



Br(B^X s l + l) 

FIG. 3: Branching ratios of B s — > fx + fi~ versus B — > X s £ + £~ in the J r -SU(5) model. 
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To see whether there are solutions simultaneously satisfied with the allowed ranges for 
these data, we can offer now some predications for B s — > £ + £~ / y, which might be measured 
at the LHCb and B factories. The numerical results are illustrated in Fig. HI We can see 
clearly that under the constraints from the inclusive decays B — > X s ^ and B — > X s £ + £~ , 
exclusive decays B s — > as well as CKM measurements extracted by the tree-level B 

decays, the branching ratio, which is very sensitive to tan (3 and charged Higgs boson mass, 
can still be up to (4 ~ 5) x 10~ 8 . Thus, it may be tested by the LHCb soon. 



CO 
CQ 



a m h * = 500 GeV 

ant, = 40 
x m„. =3000 GeV 

tanp = 2 




m u (GeV) 

FIG. 4: Branching ratio of B s — > £ + l~~{ with the combained constraints from B — > X s j, B — > 
X s £ + £~ and B s — > /i + pT . Red cross stands for the type inputs (tan/3 = 2, m h + = 3000GeV) and 
green triangle for (tan f3 = 40, m h + = 500GeV) in the F-SU (5) model, respectively. 



Rare B decays continue to be the valuable probes of physics beyond the SM. In the 
current early phase of the LHC era, the exclusive modes with muons in the final states are 
among the most promising decays. The decay B s — > fi + fi~ is likely to be confirmed before 



the end of 2012 



36|. If an enhancement beyond 10 8 and further non-decoupling effects 



are observed, we will have an indication of the J r -5't/(5) model. Although there are some 
theoretical challenges including calculation of the hadronic form factors and non-factorable 
corrections, B s — > £ + £~'-f can be expected as the next goal once B s — > fi + fi~ measurement 
is finished since the final states can be identified easily and branching ratios are large. Our 
predictions for such processes can be tested in the near future. 
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V. SUMMARY 



In this paper, we studied the vector-like quark contributions to B physics processes in 
the J r -5'[/(5) model, including the quark mass spectra, Feynman rules, the new operators 
in low energy effective theory and the correspondence Wilson coefficients, etc. As for the 
first time study, we focus on the implication of mass scale of vector like quark. The main 
conclusions we obtained are the following: 

1. There exists the sbZ interaction at tree level, and the Yukawa interactions are changed. 
The new operators 0' g and O[ must be introduced in effective Hamiltonian, and the 
Wilson coefficients are changed due to the violation of the unitarity condition. 

2. Different from many previous models, the effects of vector-like quarks on rare B decays 
such as B — > A s 7 and B — > X s £ + £~ do not decouple in some allowed parameter space, 
especially when the vector-like quark mass is comparable to the charged Higgs boson 
mass. 

3. Under the constraints from B — > A s 7 and B — > X s £ + £~, there exist scenarios in the 
model the latest measurement for B s — > fi + fi~ can be explained naturally, and the 
branching ratio of B s — > £ + £~^ can be up to (4 ~ 5) x 10~ 8 . 

All in all, due to the participation of vector-like particles, the F-SU^) model is different 
from the ordinary models such as 2HDM. In particular, the non-decouling effects are much 
more predictable and may be tested in the near future experiments. Finally, we should note 
that the large input parameter space and the sparticle effects in the F-SU^) model needs 
further work. 

Acknowledgments 

This research was supported in part by the Natural Science Foundation of China under 
grant numbers 11005006, 11172008, 10821504, 11075194, and 11135003, by the DOE grant 
DE-FG03-95-Er-40917, and by the Doctor Foundation of BJUT No. X0006015201102. 



17 



Appendix 



The loop functions for calculating the Wilson coefficients at the matching scale are the 
following 

5x + 38x 2 - 55x 2 Ax - 17x 2 + 15x 3 



A(x 
B(x 
C(x 

D(x 
P{x 
Q(x 
R(x 
S{x 
T(x 



36(x - 1) ; 
X + x 2 



X 



6(x - l) 4 
\nx, 



lnx, 



A(x-l) 2 2(x-iy 
20 - 19x + 5x 2 -2 + x 

18(x-l) 3 3(x-l) 4 
— 5 — 5x + Ax 2 2x — x 2 

12(x- l) 3 + 2(x- l) 4 n3; ' 

4(x-l) + 4(x-l) 2lnX ' 

x 2 — 6x 3x 2 + 6x 

8(x-l) + 8(x-l) 2 ' 

31x 2 + 20x 3 -4 + 18x - 30x 2 + 6x 3 



9 x- 1 



9 fx - 1 



lnx, 



38 - 79x + 47x 2 -4x + 6x 2 - 3x 4 , 

H — In X, 



108(x- l) 3 
x — 5x 2 — 2x 3 
12fx - l) 3 



18(x- l) 4 



2fx - 1] 



lnx. 
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